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ABSTRACT 


The  negative  gradient  method  is  extended  to  the  stable  analog 
computer  programming  of  a  class  of  time  varying  trajectory  problems 
defined  by  <|>(x,t)  =  0,  f(x,u,t)  =0,  u  =  x  where  x  and  u  are  n-di- 
mensional  vdctots,  (jl  is  an  m-dimensional  position  vector  function 
(m  <  n) ,  and  f  is  An  (n  -  m)  dimensional  velocity  vector  function. 
This  class  of  problems  includes  the  amplitude-stabilized  oscillator, 
two-dimensional  contour  tracing  (ref  1),  conic  section  gfeneration 
and  three-dimensional  trajectory  plotting  when  at  least  t.he  first 
integral  of  the  equations  of  motion  are  available. 


An  augmented  velocity  vector  function  is  defined 


f  '= 


f  + 


d(J> 

dt 


which  provides  n  independent  equations 
these  equations  can  be  solved  for  u  in 
the  compute  mode  u  is  connected  to  the 
resulting  system  is  analyzed  and  shown 


f  /  =  0.  For  a  fixed  x, 
the  computer  reset  mode, 
integrator  developing  x. 
to  be  easily  stabilized. 


In 

The 


1.  INTRODUCTION 

1.1  Definition  of  the  Problem 


The  problem  will  be  defined  in  n-dimensional  space.  Vector/ 
matrix  notation  will  be  employed  to  reduce  the  labor  of  writing  equa¬ 
tions.  The  following  definitions  will  be  used: 


Cl)  position  vector  x  s  [x  ,  x  ,  ...  x  ] 

1  Z  XI 

(2)  i~th  position  function  <|>.  s  (^.(x,t);  i  =  l...m,  m  <n 

(t  is  theThdependent  variable  time) 

(3)  position  function  vector  ♦  *  fV  •  •• 

(4)  velocity  vector  u  s  [u  ,  u  ,  . . „u  ] 

l  a  n 

(5)  j~th  velocity  function  f.  s  f.(u,x,t);  j  =  m  +  1,  m  +  2,  . 

J  J 


(6) 

velocity 

function  vector  f  a  [: 

£m+l ’  f m+2 ) 

(7) 

position 

error  vector 

vs  [vx, 

V*** 

(8) 

velocity 

error  vector 

a  s  f>i> 

V’" 

(9) 

gradient 

vector  V  = 
x  ! 

r  a  a 

Lax^  3x2' 

3  ’ 

1 

(10) 

gradient 

vector  V  = 
u 

[  d  ■  9 

Uu^  dug-’ 

s 

1 

(11)  general  matrix  notation  A  s  (a  )  and  I  =  unity  matrix  s 

ij 

(12)  diagonal  gain  matrix  K  s  (k^  6  )  i,j  =  1. . „n 

(13)  diagonal  gain  matrix  G  =  (g^  6^  )  i,j  =  1...  n 


.  n 


) 
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(14)  scalar  or  inner  product  of  two  n-dimensional  vectors  (p,q): 

<P,q>  =  +  p2q2  +.  ■  ■  pnqn 

(15)  positive  definite  matrix  s  a  real  symmetric  matrix  all  of 

whose  characteristic  roots  are  positive 

Using  the  aboye  definitions  the  problem  can  be  stated  as  fol¬ 
lows  :  Given  the  set  of  independent  equations 


<5>(x,t)  =  0 

(1,D 

f(x,u,t)  =  0 

(1.2) 

(1) 

.  _  dx 

U  =  X  S  dt 

(1.3) 

program  an  analog  computer  to  generate  the  unknowns  x(t),  u(t).  (1.1) 

represents  m  equations,  (1.2)  represents  n-m  equations  and  (1.3)  rep¬ 
resents  n  equations  for  a  total  of  2n  equations  in  2n  unknowns  x^,.:; 
x„,...  x  ,  u^,  u2,...  u„.  The  problem  has  been  formulated  in  terms 
of  u  and  x  rather  than  3c  and  x  because  of  the  distinction  between  the 
computer  variable,  x,  which  by  definition  is  the  total  input  to  the 
x  integrator  and  the  problem  variable  u  (i.e.  the  x  required  as  a 
solution  to  the  mathematical  problem.).  A  full  statement  of  the 
problem  also  requires  that  n-m  of  the  position  coordinates  be  speci¬ 
fied  as  initial  conditions.  Before  continuing  with  the  analysis  of 
such  a  system,  it  may  be  helpful  to  expand  upon  the  terse  problem 
statement  by  considering  some  examples  and  applications  in  two-  and 
three-dimensional  space. 

1.2  Examples  and  Applications 


In  this  section  vector  notation  is  momentarily  abandoned 
in  favor  of  more  conventional  x,y,z  notation.  Also  u  will  be  elim¬ 
inated  from  the  equations  by  the  substitution  x  =  u. 

'  Consider  the  three-dimensional  system  (2): 

(J  (x,y,z)  s  xs  +  y*  -  csz3  =  0 

<|>2 (x,y,z)  =  Cfc  +  px  +  +  ,h  0  (2) 

f3(x,y,z,x,y,z)  s  x3  +  y3  +  z3  -  s2  =  0 

(c,  a,  (3  Y,  h,  s  are  arbitrary  constants) 

The  sim,uJ.taneous  solution  of  (2)  is  a  conic  section  traced  at  constant 
speed  s'.  In  (2)  the  variable  t  does  not  appear  explicitly  in or  f. 

Consider  the  two-dimensional  system  (3): 
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(3) 


<j*^(x ,y,t)  9  x8  +  ya  -  ^(t)  =  0 

f2(x,y,x,y,t)  s  x3  +  y3  -  sa(t)  =  0 

If  r(t)  is  constant  and  s(t)  is  a  linear  ramp,  (3)  represents  linear 
frequency  modulation  without  any  amplitude  modulation.  If  both  r(t) 
and  s(t)  are  constant,  (3)  represents  an  oscillator  with  highly 
stable  amplitude  and  frequency  characteristics.  If  r(t)  is  constant 

and  s(t)  =  4r»  (3)  performs  trigonometric  resolution  (ref  2).  If  , 

QT  ' 

s(t)  is  a  constant,  (3)  represents  a  complex  modulation  scheme  obey¬ 
ing  the  law  (AM)  x  (FM)  =  constant.  The  modulating  intelligence  r(t) 
can  be  received  by  either  an  AM  or  FM  receiver  and  is  redundant. 

Consider  the  system  (4): 

(^(x^t)  s  xa  +  ya  -  rs(t)  *  0 
f2(x,y,x,y,t)  =  xa  +  y3  -  r3(t)  sa(t)  =  0 

If  in  (4)  both  r(t)  and  s(t)  vary  independently,  the  system  will  rep¬ 
resent  a  simultaneous  AM  and  FM  waveform  with  independent  messages 
r(t)  and  s(t)  . 

2.  ANALYSIS 

2.1  Partitioning  the  Problem 

Consider  first  the  situation  in  which  the  analog  computer 
is  in  the  hold  or  reset  mode.  In  this  mode  the  correct  stationary 
values  of  x  and  u  must  be  generated.  Since  <j>  is  not  a  function  of  u, 
it  is  possible  to  generate  x  by  programming  the  equations 

4x(x)  =  0,  <}2(x)  =  0,...  ^(x)  =  0  (5) 


with  n-m  of  the  coordinates  of  x  specified  as  initial  conditions. 
This,  of  course,  results  in  a  system  of  m  independent  equations  with 
m  unknowns.  Since  these  equations,  in  general,  are  nonlinear,  the 
gradient  method*  of  programming  is  required  to  guarantee  stability. 

Turning  now  to  the  generation  of  u,  the  condition  f  =  0 
provides  only  n-m  equations  to  be  solved  for  n  unknowns.  To  obtain 
the  m  additional  equations  an  augmented  velocity  function  vector  (f  ) 
must  be  defined. 


*  The  gradient  method  is  synonymous  with  least  squares,  steepest 
descent,  or  transpose  matrix  method. 


In  vector  notation 


*'  -  5t  +  1  <6> 

For  the  rigorous  interpretation  of  (6),  the  original  definitions  of 
f  and  must  be  expanded  to  n-dimensions,  with  the  first  m  components 
of  f  being  identically  zero,  and  the  last  n-m  components  of  (j>  being 
identically  zero.  Written  out: 

_  d42  ^  i 

f  =  [dt  ’  dt  dt  3  fm+l>  fm+2,,0°  fn  (  J 

If  in  the  first  m  coordinates  of  f'  the  substitution  of  u  for  x  is 
made,  then 


f'  =  0  (8) 

represents  a  system  of  n  equations  in  n  unknowns  u  ,  u  , .  .  „u  „ 

The  system  (8)  can  be  programmed  by  the  gradient  method  and  requires 
that  x  developed  from  (5)  be  inserted  as  a  parameter.  The  conditions 
obtained  in  reset  or  hold  then  will  be 

<j>  =  0 


(9)  represents  the  partitioned  system. 
2.2  Closing  the  Switch 


Refer  now  to  figure  1,  which  is  a  simplified  schematic  il¬ 
lustrating  the  connections  for  the  i-th  component  of  x  and  u.  In 
reset  or  hold,  switch  S  is  open  and  u  and  x  assume  their  correct 
stationary  values.  It  is  reasonable  to  suppose  that  if  S  were  closed, 
the  resulting  system  would  produce  a  very  close  approximation  to  the 
desired  trajectory,  provided  the  system  maintains  dynamic  stability. 
The  resulting  computer  differential  equation  will  now  be  derived. 


In  accordance  with  the  gradient  method. 


let 

V  = 

-  grad^ 

«M>  = 

-vx<M> 

(10) 

a  = 

-  grad^ 

(f',f ') 

=  -Vu(f',f') 

(ID 

From  figure  1 

we  have 

x  = 

OH  +  Ky 

(12) 

u  = 

Ga 

(13) 
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(a  is  the  potentiometer  setting,  figure  1,  and  is  the  time  scale 
factor)  Differentiating  (12)  with  respect  to  time  and  substi¬ 
tuting  (13) 

x  =  QGa  +  Kv  (14) 


The  system  stability  may  be  investigated  by  linearizing  the  system 
about  some  arbitrary  operating  point .and  examining  the  effect  of 
small  perturbations  6x  and  §u.  Operating  on  equations  (12)  and 
(14)  with  the  variational  operator  6 


d(6x) 

dt 

d2(6x) 

dt* 


as5u  +  K6v 

:  OG6a  +  K 


d(6v) 

dt 


(15) 

(16) 


It  is  now  necessary  to  express  6v  and  6a  in  terms  of  6x.  Taking 
variations  of  both  sides  of  (10) 


(17) 


Since  variations  are  taken  starting  from  an  assumed  equilibrium 
state  ^  =  0  the  second  term  on  the  right  in  (17)  vanishes.  The 
matrix 


is  positive  definite  since  it  can  be  factored  into  the  product  of 
a  matrix  by  its  transpose;  (17)  then  becomes 

6v  =  -A  6x  (18) 

Note  that  since  <j)  is  not  a  function  of  u,  only  the  variation  with 
respect  to  x  had  to  be  considered  in  (17).  To  .'find  6a>  take  varia¬ 
tions  of  both:  sides  of  (11),  noting  now  that  both  gu  and  6x  will 
contribute  to  §a. 
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Figure  1.  Simplified  computer  diagram  showing  connections  for 


In  (19)  terms  involving  second  partial  derivatives  are  absent  be¬ 
cause  of  the  assumption  that  the  initial  conditions  <)>  =  0  and  f '  =  0 
are  satisfied. 


Rewriting  (19) 


6a  =  -B  6u  -H  6x 


whe  re 


B*  <V 


H=  <V 


(20) 


The  matrix  B  is  positive  definite  but  the  matrix  H  is  not  so  re¬ 
stricted.  Substituting  (18)  and  (20)  in  (16) 

j  2  j 

^2  (6x)  +  aGBfiu  +  aGHfix  +  K  —  (A6x)  =  0  (21) 

d3  riA  d 

^  (6x)  +  aGB6u  +  0GH6x  +  K  ^  6x  +  KA  —  (6x)  =  0  (22) 

From  (15)  and  (18) 

Os5u  =  KA6x  +  ~  (fix)  (23) 

dt 

Substituting  (23)  in  (22) 

ha 

I  -zzs  (6X>  +  (KA  +  GB)  -r£  (fix)  +  (GBKA  +  OSGH  +  K  ~)  fix  =  0  (24) 

at  at  ax 

(24)  is  the  characteristic  matrix  differential  equation  of  the  sys¬ 
tem.  Note  that  (I)  the  unity  matrix  is  positive  definite  and 
(KA  +  GB)  is  positive  definite ,  since  it  is  the  sum  of  two  positive 
definite  matrices  KA  and  GB.  Now  it  can  be  shown  by  an  extension 
of  the  argument  of  Bellman  (ref  3)  that  (24)  will  be  stable  if  the 
matrix  sum 

dA 

GBKA  +  OGH  +  K  —  (25) 

dt 

can  be  expressed  as  the  sum  of  a  positive  definite  matrix  and  a  skew 
symmetric  matrix.  The  matrix  product  GBKA  is  positive  definite  and 
dA 

and  the  matrix  K  —  is  symmetric „  Let  H  be  expressed  as  the  sum  of 
dt 

a  symmetric  part  H*  and  a  skew  symmetric  part  S 
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H  =  ff  +  S 

Let  the  notation  A  >  B 
that  A  -  B  is  positive 
can  be  written 

GBKA  +  qGH  + 

or 

GBKA  >  -  QGH 

It  is  now  evident  that 
scale  factor  CL  must  be 
than  the  real  time  case  (a  =  1).  If  the  term  K  ^  is  causing  in" 
stability  K  can  be  decreased.  If  only  the  full  set  of  velocity 
functions  is  given  (f  ^  f^,,,.  f  ),  then  A  =  0  and  the  condition 
for  stability  becomes  JT  >  0.  TBis  may  be  impossible  to  satisfy, 
in  which  case  at  least  one  function  (j> .  must  be  found  by  integration 
from  the  system  f  =  0.  Since  at  lea^t  (^.  is  available,  A  is  rein¬ 
stated  and  System  stability  is  obtained  a1^  before. 

3.  CONCLUSION 

It  has  been  shown  that  the  powerful  negative  gradient  technique 
can  be  extended  to  successfully  program  trajectory  problems  as  de¬ 
fined  in  the  introduction.  Computer  stabilization  may  require  time 
scaling  (decreasing  a)  or  reducing  the  gain  K. 
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>0 


(27) 


dA 
K  dt 


(28) 
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